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Abstract. The problem of determining the pattern of gas motions in the central regions of 
disk spiral galaxies is considered. Two fundamentally different cases - noncircular motions in the 
triaxial bar potential and motions in circular orbits but with orientation parameters different from 
those of the main disk ~ are shown to have similar observational manifestations in the linc-of-sight 
velocity field of the gas. A reliable criterion is needed for the observational data to be properly 
interpreted. To find such a criterion, we analyze two-dimensional nonlinear hydrodynamics mod- 
els of gas motions in barred disk galaxies. The gas line-of-sight velocity and surface brightness 
distributions in the plane of the sky are constructed for various inclinations of the galactic plane 
to the line-of-sight and bar orientation angles. We show that using models of circular motions 
for inclinations i > 60° to analyze the velocity field can lead to the erroneous conclusions of a 
"tilted (polar) disk" at the circumnuclear region. However, it is possible to distinguish bars from 
tilted disks by comparing the mutual orientations of the photometric and dynamical axes. As an 
example, we consider the velocity field of the ionized gas in the galaxy NGC 972. 
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1. Introduction 

The analyses of the structure and kinematics of 
spiral galaxies often require detailed information 
on the gas motion patterns in the stellar systems 
considered: the rotation curve and whether the 
velocity field includes noncircular motions. If gas 
motions are circular it is enough to measure line- 
of-sight velocity profiles along two sections pass- 
ing through the disk center to unambiguously 
infer the orientation parameters and determine 
the rotation curve (Zasov 1993). This however 
only a rough approximation of the real dynamics 
of the gaseous disk: the velocity fields of spiral 
galaxies can exhibit large-scale noncircular mo- 
tions amounting to 20 — 30kms~^ , which are 
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due by the spiral pattern. 

Much more detailed information can be ob- 
tained by measuring the full velocity field, i.e., 
by determining the line-of-sight velocity distribu- 
tion over the entire galaxy disk. This has since 
long become common practice among radio as- 
tronomers who use interferometry technique to 
measure the spectral lines of neutral or molec- 
ular gas. Similar results are obtained at optical 
wavelengths with 2D-spectroscopy using Fabry- 
Perot interferometers or integral field spectro- 
graphs. The common feature of all these meth- 
ods is that they all involve construction of the 
so-called "data cube" (Tully 1974). A special 
technique allows the observational data to be re- 
duced to a form where each image element has 
its own individual emission- (absorption-) line 
spectrum. The line-of-sight velocity field is then 
constructed from the Doppler shifts of the lines 
studied. 
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It is, however, often difficult to unambigu- 
ously reconstruct the rotation curve and ana- 
lyze the direction of noncircular motions in the 
galaxy even if the full line-of-sight velocity field 
is available. Noncircular motions, which are pri- 
marily due to the spiral pattern, make it impos- 
sible to isolate the circular rotation component 
simply by averaging the total gas velocity field. 
Fourier analysis of the line-of-sight velocity field 
(Lyakhovich et al. 1997) allows the full spatial 
pattern of gas motions in a spiral galaxy to be 
reconstructed under certain assumptions about 
the nature of the spiral structure. 

Barred galaxies exhibit noncircular motions 
of even larger amplitude. Although the interpre- 
tation of observed gas velocities in galaxy bars re- 
mains highly ambiguous, a number of important 
conclusions to be made about the motion pattern 
in the bar can be drawn by analyzing the gas ro- 
tation in terms of the circular-motion model. The 
bar produces a turn of the dynamical axis (i.e., 
the line of maximum line-of-sight velocity gradi- 
ent) relative to the line-of-nodes (Chevalier and 
Furenlid 1978; Afanasiev et al. 1992; Zasov & 
Moiseev 1999) and this turn can be determined 
using the well-known algorithms for the analysis 
of disk galaxy velocity fields. 

The most comprehensive conclusions about 
a bar in a galaxy can be drawn from numerous 
works on nonlinear computer simulations involv- 
ing hydrodynamics equations written for "gas" 
consisting of macromolccules in the form of gas 
clouds. See Lindblad (1996) for a review of var- 
ious model simulations. The common result of 
these works is that a bar to be an efficient mecha- 
nism for removing the angular momentum of the 
gas, thereby making the latter to partly move 
centerward along the bar and lose the energy 
at its shock edges (Athanassoula 1992a, 1992b; 
Levy et al. 1996). This motion pattern distorts 
the observed velocity contours, aligning them 
along the bar, whereas the dynamic axis, which 
is perpendicular to the isovelocities it crosses, 
must, accordingly, turn in the opposite 

We must, nevertheless, bear in mind that the 
turn of the dynamic axis can also be caused by a 
tilted disk at galaxy center (Zasov & Sil'chenko 



1996; Zasov & Moiseev 1998). However, in this 
case, the axis of a circular rotating disk coin- 
cides with the line-of-nodes. Also the photomet- 
ric and dynamic axes should therefore turn in 
the same direction. These are only qualitative 
considerations; therefore, concrete computations 
are needed, because the projection effect with ra- 
dial gas flows in the bar can produce a pattern 
similar to that expected in the tilted disk case. 

Extensive panoramic spectroscopy of HII ve- 
locity fields in spiral galaxies have been obtained 
with the Fabry- Perot interferometer attached to 
the 6-m BTA telescope of the Special Astro- 
physical Observatory of the Russian Academy of 
Sciences. It is very important, when interpret- 
ing these data, to analyze gas motions in the 
bars and identify eventual minibars (smaller than 
1 kpc) from the dynamical manifestations they 
cause in the velocity field (see, e.g., Afanasiev et 
al. 1989; Zasov & Sil'chenko 1996; Sil'chenko et 
al. 1997). To investigate whether the above pat- 
terns in the behavior of the mutual orientations 
of photometric and dynamical axes can in prin- 
ciple be applied to analyze observational data, it 
would be useful to consider the results of numer- 
ical velocity-field simulations with allowance for 
sky-plane projection effects and limited spatial 
and spectral resolution. 

In this paper, we use the initial data, which 
are similar to the results of computations by 
Levy et al. (1996), to construct simulated HII 
velocity fields similar to those actually observed 
in the Hq and [Nil] lines with the standard 
Fabry-Perot attached to the BTA telescope (see 
Dodonov et al. (1995) for a detailed description 
of the instrument). These data have a spatial 
resolution of (1.5 — 4)" and a spectral resolu- 
tion (50 — 150) km s^^ in velocity terms, with 
an accuracy of inferred line-of-sight velocities of 
(3- 10)kms-i . 

2. Velocity fields constructing 

Our simulations are based on the numerical so- 
lution of hydrodynamics equations in an exter- 
nal gravitational field. The adopted gravitational 
potential model consists of a nonaxisymmetric 
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perturbation superimposed on an axisymmetri- 
cal component (bulge, disk, and halo) see Mat- 
suda et al.(1987). 

The nonaxisymmetric perturbation (stellar 
bar mode) is turned on smoothly and gradually 
increased until it reaches the fixed level. The 
model exhibits a specific quasi-periodic regime 
of gas passage through the bar (Levy et al. 1994, 
1996). Interestingly, the resulting flow pattern is 
bisymmetricthe gas flows centerward along the 
bar in two sectors, whereas in the other two sec- 
tors the flow direction is reversed; i.e., the gas 
flows away from the disk center due to the saddle 
point of the gravitational potential. The develop- 
ing ar also produces a two-armed spiral pattern 
outside the corotation circle of the bar. 

We used the results of simulations identi- 
cal to those described in detail by Levy et al. 
(1996). The only difference was that we party 
adjusted the gravitational potential to closely re- 
produce the rotation curve of our Galaxy (Haud, 
1979). Our model had a maximum linear rota- 
tion velocity of 250kms~^ ; disk scale length 
of 3 kpc; radial bar scale length of 1 kpc, and 
a radius of a simulated region 10 kpc. We fur- 
ther assumed that the galaxy is at a distance 
of 20.6 Mpc, where 1" corresponds to 0.1 kpc, 
and has a systemic line-of-sight velocity (i.e., the 
center-of-mass velocity relative to the observer) 
of Vsys = 1545 km s~^ . 

Our numerical simulations yielded the distri- 
butions of gas surface density ayas^R, f), radial 
VR{R,ip) and azimuthal V^(i?, 99) gas velocity 
components. Here, R and ip are the radial and az- 
imuthal coordinates in the galactic plane, respec- 
tively. The polar grid cells had sizes of A(p = 2° 
in the azimuthal and 0.016 < AR < 0.2 kpc in 
the radial direction, depending on the galacto- 
centric distance. Based on these data, we con- 
structed model cubes using the following algo- 
rithm. Line-of-sight velocity at point {R, (p) : 



Vobs (R, ^) = Vsys + Vr{R, if) sin if sin i 

+ V^{R, if) cos If sin i 



(1) 



where Vgys and i arc the systemic linc-of-sight 
velocity and the inclination of the galaxy plane 
to the sky plane, respectively, ip = set for the 



line-of-nodcs. 

We projected distributions VobsiR,v) 
a gas {R, f) onto a Cartesian grid with a cell size of 
0.02 kpc (O'.'S), assuming that the position angle 
of the line of nodes in the sky plane is PAq = 90° . 
As a result, we obtained the distributions of 
line-of-sight velocity and gas surface density dis- 
tributions Vohs{x-,y) and agasix,y), respectively. 
The emission-line spectrum at each point is flt- 
ted by a Gaussian centered on Vobsix,y) with 
a FWHM = 130 km s"^, which approximately 
corresponds to the spectral resolution of the 
Fabry-Perot Interferometer attached to the BTA 
telescope. We neglected the velocity dispersion 
of individual gas clouds during a construction of 
the total spectrum because of the coarse spectral 
resolution. We assumed that the spectral-line in- 
tensity is proportional to (Jgasix^y). The images 
in all spectral channels were Gaussian smoothed 
to make the resulting spatial resolution equal to 
2", which is close to the typical seeing value dur- 
ing real observations. 

Based on the smoothed-cube spectra, we con- 
structed the velocity field and the spectral-line 
image (the observed distribution of gas surface 
density) and determined the observed orienta- 
tion parameters of the gaseous disk using the 
technique similar to that described by Begeman 
(1989). A similar algorithm is used in the well- 
known GIPSI radio astronomical data reduction 
software to analyze the observed velocity fields. 

Below, we briefly describe the procedure used 
for the analysis. 

Introduce polar coordinate system (r, PA) 
in the sky plane, where r and PA are the dis- 
tance from the rotation center and position an- 
gle, respectively. The observed line-of-sight ve- 
locity Vobs ill case of purely circular rotation 
is 

Vobs{r,PA)=Vsys 

A-V JFKr)) cos(PA-PAo)sint [2) 

-r Vrot{-n{T)) (i+si„2(p^_p^^)t^n2 j)i/2 

and the distance from the rotation center in 
the galaxy plane is 

R{r) = r(l + sm^{PA - PAq) tan^ z) V2 (3) 
where Vrat and PAq are the circular rotation 
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velocity and position angle of the line-of-nodes, 
respectively. 

During reduction, the observed velocity field 
is divided into elliptical rings defined by equa- 
tion (3) for R = const. A nonlinear least squares 
technique are used to fit model curve (2) to the 
observed dependence Vobs{P-^)- As a result, we 
obtain for each radius r the corresponding disk 
orientation parameters PA and i, and the circu- 
lar rotation velocity Vrot- 

Note that while the positions of the ex- 
tremes of function (2) allow PA to be determined 
quite unambiguously; i estimates are much more 
uncertain. Thus, at small inclination angles i 
(galaxies seen as face-on), it is impossible to 
separate the contributions of Vrot and i to the 
observed velocity projection and only Vrot sin i 
can be unambiguously inferred. Begeman (1989) 
used a reduction of radio data to show the cou- 
pling of the two parameters at inclination angles 
i < 40°. 

We constructed and analyzed simulated fields 
for various bar development stages and inclina- 
tion angles ranging from 30° to 70°. We varied 
the angle between the linc-of-nodes and the bar 
with a step of 15° and constructed a total of 
more than 400 simulated galaxies. 

3. Analysis of the results 

In Fig. 1, we plot the gas surface-density dis- 
tribution and line-of-sight velocity contours for 
various bar orientation angles. This is a typical 
pattern for dimensionless time instants T < 3.4 
(where T = 1 corresponds to one bar revolution) . 
The bar then saturates and the spirals degener- 
ate into a pseudo-ring (see Levy et al. (1996) for 
details). The contour turns inside r < (10 — 15)" 
are caused by the bar, but on greater galacto- 
centric distances they are determined by the mo- 
tions in spiral arms. The typical S-shapcd radial- 
velocity contours can be seen in the central re- 
gion. It is interesting that the contours in ques- 
tion are distorted even if the bar is aligned along 
the major axis of the galaxy (0 = 0, where (j) is 
the angle between the major axis of the bar and 
the line-of-nodes measured in the galaxy plane). 



The line-of-sight velocity projection is formally 
equal to zero in this case; however, the bar has 
a finite width and velocity perturbations contain 
an azimuthal component, which also contributes 
to the radial velocity. 

Fig. 2 shows typical results for dynamical axis 
orientation for an inclination angle of i = 70°. 
At the very center of the galaxy, the dynamical 
axis always deviates significantly from the line of 
nodes PAq = 90° (except ^ «i 0°, -90° i.e., in 
the vicinity of the major and minor axes, respec- 
tively). Beginning with r < 5", the angle between 
the two lines decreases and becomes negligible 
at r « 10" (in the vicinity of the bar tips). At 
greater galactocentric distances, the dynamical 
axis turns to the other side of the line-of-nodes 
and misalignment disappears only at r > 20". 

As the bar evolves, a shock develops at its 
leading edge and this shock is especially con- 
spicuous in the pressure maps due to the strong 
gas density in the region considered. Such an 
arrangement of shock fronts is consistent with 
both the observations of dust lanes in bars and 
with the results of earlier numerical simula- 
tions (Athanassoula 1992a, 1992b). The beat- 
frequency modulation about PAq is due to the 
complex velocity pattern at the bar tips, where 
the shock front moves from the leading bar edge 
to the inner edge of the trailing arms (because 
the corotation radius is r f« 20"). The PA vari- 
ation pattern is asymmetric with respect to ^ = 
—90°, because of a preferred direction of angle- 
counting determined by the pitch angle of the 
spiral density wave, whose direction is opposite 
to that of galactic rotation (trailing spiral). 

The thick line in Fig. 2 shows the orienta- 
tion of the sky-plane projection of the (the ob- 
served photometric) axis of the bar. It is evident 
from the figure that, in spite of the complex pat- 
tern of PA radial variation, the dynamical axis 
in the bar (r < 10") always turns with respect 
to the photometric axis in the direction opposite 
that of the line-of-nodes: the two lines turn in an- 
tiphase. This conclusion applies to all simulated 
fields with different i, (j) and T. 

In order to better illustrate this effect. Fig. 
3 shown the deviations of the dynamical axis 
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Figure 1: Sky plane projection of the gas surface brightness in the central region for various bar turn 
angles (p (indicated at the top right corner in each graph) at time T = 2.0 since the start of the 
simulation. Contours show the Vgys and Vgys ± 100 km s~^ line-of-sight isolines of model velocity field. 
The galaxy plane tilt is i = 70°; the dashed curve shows the line-of-nodes (position angle PA = 90°). 



(averaged over all points with r < 6" depends 
on the bar orientation. In this figure, the am- 
plitudes of the turn of the dynamical axis are 
somewhat smoothed by averaging; however, the 

overall trend remains conspicuous at all inclina- 
tion angles and various bar development stages. 

At angles ranging from 0° to ±45° (which 
correspond to PA pa ±(30° — 45°) depending 
on the inclination angle i), the dynamical axis 
turns virtually in antiphase with the bar orien- 
tation. With a further increase of the angle be- 
tween the bar and the line-of-nodes the dynam- 
ical axis begins to "lag behind"; its orientation 
changes slower than that of the photometric axis 
and it virtually coincides with the linc-of-nodcs 
at ~ ±90°. However, the turn of the dynamical 
axis from the line-of-nodes is always in the sense 



opposite that of the bar axis. The lag of the dy- 
namical axis is a result of a combined effect of 
projection and that the bar fails to stop the gas 
rotation completely (the maximum bar contrast 
in our experiment was 30%), and there is always 
an azimuthal component besides the radial mo- 
tions. 

We now analyze the behavior of the inclina- 
tion i formally inferred from the velocity fields 
(Fig. 2). Radial motions in the bar region distort 
relation (2), which should have been observed if 
the rotation is purely circular. Therefore the in- 
clination inferred at these radii can differ sub- 
stantially from the true inclination of the disk. 
The variations of i can have amplitudes as high 
as 20° — 30° and are often rather chaotic. How- 
ever, the following pattern can be observed: for 
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Figure 2: Results for model velocity fields corresponding to various bar turn angles at time T = 2. 
Radial variations of the dynamical axis orientation PA and tilt angle i. Solid lines in PA plots show 
the position angle of the bar projected on the sky plane (indicated by "PAbar = on each graph). 
Dashed lines on i plots show the galaxy inclination i = 70° 



galaxy inclinations i > 60° — 70°, the formally 
inferred inclination in the bar region becomes 
virtually constant at ^ > -30° {(p < -150°). 
One gets a false conclusion that the galaxy cen- 
ter contains a decoupled disk whose inclination 



to the sky plane is 15° — 25° smaller than that of 
the " outer" disk of the galaxy. The real angle be- 
tween the disk planes can be as high as 90° (the 
polar disk). 

The point is that, although at i > 50° - 60°, 
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Figure 3: Mean deviation of the position angle of the hne-of-nodes in the bar region APA^iyn plotted 
as a function of the bar position angle APAf,ar for i = 70°, 60°, 30°. Different signs correspond to 
different dimensionless time instants T. 



the Vobs dependence on PA in the equation (2) 
differs strongly from the cosine relation, radial 
motions whose projected velocities are maximum 
at the minor axis of the galaxy distort the pat- 
tern implied by equation (2) making it more like 
a cosinusoidal corresponding to smaller i. It is 
thus evident that restricting the analysis to sur- 
face photometry alone makes it rather difficult to 
distinguish a minibar from a nuclear polar disk; 
one may erroneously take the observed pattern 
for a manifestation of gas moving in curricular 
orbits with smaller inclination angles, whereas it 
might simply be a line-of-sight projection of a su- 
perposition of circular and noncircular motions 
in the bar. A turn of the dynamical axis should 
also be observed in the case of a real tilted (in- 
cluding polar) disk, and analyzing the velocity 
field alone can lead to wrong conclusions. Recall 
that in the tilted disk case, both the dynamical 
and photometric axes should turn in the same 
direction (Zasov & Moiseev 1999). A tilted disk 



can therefore be distinguished from a bar only 
by comparing the PA inferred from the velocity 
field analysis with the surface photometry data. 

4. Kinematics of gas in the galcLxy 
NGC 972 

In order to illustrate our approach, we applied 
it to the NGC 972 galaxy. The galaxy are ob- 
served in the [NIIJA6583 emission line with the 
scanning Fabry- Perot Interferometer attached to 
the 6-m BTA telescope of the Special Astrophys- 
ical Observatory of the Russian Academy of Sci- 
ences. The seeing and spectral resolution were 
l'.'5 and '-^ 50kms~^, respectively. Two veloc- 
ity fields with different image orientations rel- 
ative to the detector (two-dimensional photon 
counter IPCS) are constructed. The results of 
the velocity-field analysis arc briefly described by 
Zasov &; Moiseev (1999), and a more detailed pa- 
per is now in preparation. Fig. 4 shows the radial 
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profiles of the position angle of the dynamical 
axis PA and inclination angle i formally inferred 
from the velocity field in terms of the circular ro- 
tation approximation as described in Section 2. 
Dots with error bars are mean values averaged 
over both velocity fields. 

The mean inclination of the galaxy to the 
line of sight is i ~ 64°; the position angle 
of the line-of-nodes, PA = 150°, in agreement 
with the orientation of the outermost K-band 
isophotes (2.2 /x) based on observations taken 
at the UKIRT infrared telescope (Hawaii). The 
abrupt change of PA and i at r < 10" can be in- 
terpreted in terms of both inclined disk or mini- 
bar models. 

However, comparing the turns of the dynamic 
axis with the K-band isophote orientation allows 
this ambiguity to be resolved. The black line seg- 
ment in Fig. 4 indicates the mean orientation in 
the central isophotes of the galaxy, i.e., the po- 
sition of the photometric axis. It is evident from 
the figure that the orientation of the photometric 
axis (relative to the line-of-nodes) varies in an- 
tiphase with that of the dynamical axis, a pattern 
that is indicative of a ~ 10" (1 kpc) bar whose 
axis makes an angle of ^ —120° to the line-of- 
nodes. Variations of angle i in this case are due 
to noncircular motions in the bar, as noted above 
in Section 3. 

We constructed model velocity fields with ori- 
entation parameters equal to those of the NGC 
972 disk for various time instants T correspond- 
ing to a gradual increase of the contrast of the 
gravitational potential of the bar (Levy et al. 
1996). Solid lines in Fig. 4 show the results for 
the model velocity field at T = 2.2. Note that 
the observed plots agree well with the computed 
ones given the fact that, strictly speaking, the 
mass distribution in NGC 972 differs from that 
of the adopted model. 

5. CONCLUSION 

Our analysis of model velocity fields constructed 
from numerical nonlinear two-dimensional sim- 
ulations of gas flows in bars showed that, if 
the galaxy plane is inclined by more than i > 
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Figure 4: Results of an analysis of the [Nil] 
line velocity fields of the NGC 972 galaxy (black 
dots with error bars). Two large circles show the 
results of Zasov &; Silchenko (1996). Above: ra- 
dial variations of the dynamical axis orientation, 
thick interval shows the mean orientation of cen- 
tral K-band isophotes of NGC 972 (bar position). 
Below: gaseous disk tilts as inferred from veloc- 
ity fields. Solid lines in each graph show model 
dependencies for the dimensionless time instant 
T = 2.2 and bar turn angle (j) = 120°. 

60° — 70° to the sky plane, the observed distri- 
bution of gas line-of-sight velocities in the bar 
region can create an illusion of a "tilted nuclear 
disk" , implying a wrong interpretation of the ob- 
served velocity field. 

However, in the case of a bar, the dynami- 
cal axis always turns with respect to the line- 
of-nodes in antiphase with the photometric axis, 
thereby allowing the bar to be distinguished from 
a tilted disk by a combined analysis of the kine- 
matics and photometry. We illustrated the use 
of this criterion by applying it to the NGC 972 
galaxy 
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